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ABSTRACT 

The detection of Pop III supernovae could directly probe the primordial IMF for the first time, 
unveiling the properties of the first galaxies, early chemical enrichment and reionization, and the 
seeds of supermassive black holes. Growing evidence that some Pop III stars were less massive than 
100 Mq may complicate prospects for their detection, because even though they would have been more 
plentiful they would have died as core-collapse supernovae, with far less luminosity than pair-instability 
explosions. This picture greatly improves if the SN shock collides with a dense circumstellar shell 
ejected during a prior violent LBV type eruption. Such collisions can turn even dim SNe into extremely 
bright ones whose luminosities can rival those of pair-instability SNe. We present simulations of Pop 
III Type Iln SN light curves and spectra performed with the Los Alamos RAGE and SPECTRUM 
codes. Taking into account Lyman-alpha absorption in the early universe and cosmological redshifting, 
we find that 40 M Q Pop III Type Iln SNe will be visible out to z ~ 20 with JWST and out to z ~ 7 
with WFIRST. Thus, even low mass Pop III SNe can be used to probe the primeval universe. 
Subject headings: early universe - galaxies: high-redshift - stars: early-type - supernovae: general - 
radiative transfer - hydrodynamics - shocks 



1. INTRODUCTION 

The first stars in the universe are thought to form in 
10 5 - 10 6 M Q cosmological halos at z ~ 20 - 30. Unfortu- 
nately, because these stars lie at the edge of the observ- 
able universe there are no observational constraints on 
their properties. The original numerical simulations of 
primordial, or Pop III, star formation suggest that they 
are very massive, 100 - 500 Mq and that they form in 
isolation, one per halo (IBromm et al.l [19991: lAbel et al.| 
2000L 120021: IBromm et all 120021: iNakamura fe Umemura 
20011 IQ'Shea fe Norman] 120071: IWise fc Abell I2007t 
Q'Shea k Normanl 120081: lYoshida et al.l 120080 . Newer 
calculations have since found that som e Pop III stars 
may form in binaries (jTurk et al.l 120091) or even small 



swarms of 20 - 40 Mq stars dStacv et al.l 



2010; 



2011 



Clark et al.l 12011 ISmith et al.1 120111: IGreif eTaf] 
20121) . Simulations of UV breakout from primordial 
star-forming disks suggest that ionizing feedback in 
some cases may limit the masses of the first stars to 
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20 - 50 Mp dHosokawaet al.l 120111: iStacv et al.l I2012at 
Hosokawa et al.l 120121) fb ut see also | Omukai fe Pallal 



20011: IQmukai k Inutsukal l2Tj02t IcWkai fe Fallal I2003t 
Tan fc McKed 12004 IMcKee fc Tanll2008D . However, all 



these estimates must be regarded to be very preliminary, 
since no simulation has evolved a newly formed Pop III 
protostellar disk all the way to the end of the life of one 
of its stars with realistic p hysics (for recent reviews, see 
lGlovenl2rjT2] iWhalerJ 120 12ft Furthermore, the impact of 
turbulence (jLatif et al.ll2013h . magnetic fields due t o the 
small scale turbulent dynamo (|Schober et al.l 12012ft . and 
radiation transport on the evolution and stability of the 
disk are not well understood. 

There have been attempts to constrain the Pop III 
IMF by modeling the nucleosynthctic imprint of primor- 
dial supernovae (SNe) on later generations of stars. This 
imprint is now sought in the fossil abundance record, 
the pattern of chemical elements found in ancient, dim 
metal-poor stars now being surveyed in the Galactic 
halo fe.g.. iCavrel et al.l [2001 iBeers fc Christliebl 120051 
iFrebel et al.l 120051 : lLai et al.l 120081 : iCaffau et al.l 12012ft . 
Recent simulations indicate that 15-40 Mq Pop III core- 
collapse (CC) SNe may have contributed significantly 
to early chemical enrichment, further cor roborating the 
existe nce of lower-mass Pop III stars (|Joggerst et al.l 
12010ft . Some have taken the absence of the distinctive 
odd-e ven nucleosynthetic patt ern of pair-instability (PI) 
SNe (jHeger fc Wooslevl 120021 ) in extremely metal-poor 
stars to imply that there were no very massive Pop 
III stars. However, evidence of the odd-even effect has 
now been found in high-red shift damped Lyman alpha 
absorbers (Cooke e t al.ll20ill) and perhaps in a new sam- 
ple of stars from the Sloan Digital Sky Survey (SDSS, 
iRen et al~ll2012f ). It is also now known that Pop III PI 
SNe could easily have enriched later stars to metallicities 
above those targeted by surveys to date dKarlsson et al.l 
l2lMIJoggerst fc Whaler] MHi IChen" et alJl2011ft . Much 
remains to be understood about how metals from the 
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first stars are take n up into later generations by cosmo- 
logical flows (e.g. JGreiLetaD |200a IChiaki et al.l I20T1 
iRitter et al.l 120121 ) . Nevertheless, it is clear that low- 
mass Pop III stars may have been common in the early 
universe, with profoun d consequences for t h e character 
of pr i mitive galaxies (iJohnson et al.l 20081: IGreif et"aT 



2008) I Johnson et all 120091: IGreif et al.l 12010: Ueqn et al 
IT 



2011 IPawlik et all 120111 L'OUt l\Vi.-e et all 12012ft. early 



rcionization (IWhalen et all 120041: IKitavama k Yoshidal 



2T)0l t I Alvarez et alJl200fit lAbel et al. || 2007t IWise k Abell 
2008 ) and chemica l enri c hment (|Mackev et al.l 12003c 



Smith k Sigurdssonl I2M7I : iSmith et al. l 120091) . and the 



origin s of supermassive bl ac k hol es jBromm k Loeb 
2003t Uohnson k Bromml 120071: iDiorgovski et al 



2008t iMilosavlievic et all 120091: lAlvarez et al.l 120091: 



Lippai et all 120091 ITanaka fe Haimanl 120091: ILi 
20111: iPark k Ricottil 120111 l2012ar Uohnson et al 



2012al IWhalen k Frverl 120121: lAgarwal et all 120121 
Johnson et al. H2012bt IPark k Ricottill2012bl) . 



The best prospects for determining the masses of 
Pop III stars in the near term lie in detecting their 
SNe (iBromm et all [ 20031: IKitavama fc Yoshidal 12005c 
IWhalen et al.ll2008d : IVasiliev et al.ll2012l). Ev en though 
they are extremely luminous (jSchaereiTl2002[ ). individ- 
ual primordial stars arc still too dim to b e found by the 
Jame s Webb Space Telescope (JWST) (Gardner ct al. 
120061 ) or 30-meter class telescopes (although in princi- 
ple their H II regi ons could be detected via strong grav- 
itational lensing; iRvdberg et~aTl 120131 ). Pop III SNe 
can be 100,000 times brighter than their progenitors 
or the host galaxies in which they reside, and their 
masses can be inferred from their luminosity profiles. 
The main obstacle to their detection is Lyman absorp- 
tion by neutral hydrogen prior to the era of reion- 
ization, which absorbs or scatters most photons from 
these ancient explosions out of our line of sight. Pop 
III SNe must emit enough luminosity below the Ly- 
man l imit to be obs e rved in t he ne ar infrared (NIR) lo- 
cally. IWhalen et ail (|2012allbl . I2013D recently found that 
JWST will detect Pop III PI SNe at any epoch and the 
Wide-Field Infrared Survey Telescope ( WFIRST) and 
Wide-field Imaging Surveyor for High-Redshift ( WISH) 
will find th em out to z ~ 15 - 20 in all-sky NIR 
surveys (see iScaimapieco et alll2005t iKasen et all 120111; 
Pan et all l2012al lbl; iHummel et all 120121: iDessart et all 



20131 for past studies of PI SNe and their detection). 



Their extreme brightnesses make P I SNe ideal probes of 
the earliest stellar populations (see lGal-Yam et aDl2009t 
lYoung et al1l2010L on the detection of the SN 2007bi, a 
PI SN candidate in the local universe). 

Pop III CC SNe may be more plentiful at high red- 
shifts but they are more difficult to detect bec ause of 
their lower luminosities. IWhalen et all (|2012d ) calcu- 
late detection limits of z ~ 10 - 15 for 15-40 Mq Pop 
III CC SNe for JWST for explosion energies of 1 - 2 x 
10 51 erg, so they can be found in primitive galaxies but 
not in the first star-forming halos. Core-collapse explo- 
sions also cannot be used to differentiate between pri- 
mordial and Pop II progenitors in protogalaxies because 
their central engines are not very sensitive to met allicity 
(IChieffi k Limongill200l IWooslev fc Hegerll2007t) . Such 
events can therefore trace star formation rates in the first 
galaxies but are not ideal probes of the primordial IMF. 

However, some CC SNe, Type Iln SNe, have recently 



been discovered with luminosities that rival those of far 
more powerful explosions. SN 2006tf and SN 2006gy, 
whose bolometric luminosities exceed 10 44 erg s _1 and 
are on par with tho se computed for Pop III PI SNe (see, 
e.g.. lGal-Yam|[20Tl ). have now been observationally and 
theoretically connected to the collision of SN ejecta with 
a dense circumstellar shell ejected by a violent luminous 
blue variable (L BV) eruption a few years prior to the 
death of the star (ISmith fe McCravj|2007l; IGal-Yam et al l 



20071 IGal-Yam fc Leonard! 120091 : Ivan Marie et al.l I2TTT0F 



Chevalier fc Irwinl 120111 ). The shell is thought to be 



opaque because only photons from the collision are ob 
served, not those from shock breakout from the surface 
of the star. Type Iln SNe likely occur at high red- 
shifts because 15 - 40 Mq Pop III stars that do not 
have much convective mixing over their lifetimes are 
known to die as compact blue giant s (|Scannat>ieco et all 
20051: Uoggerst et al.ll20ioh (see also lEkstrom et al.ll200S : 



Stacv et all 120111: lYoon etaTI 12012c IStacv et al.l 12012b . 
on the effect of rotation and magnetic fields on the 
evolution of Pop III stars). They may be observable 
at redshifts above those at which normal luminosity 
CC SNe can be detected because the shell becomes ex- 
tre mely luminous in UV when the ejecta crashes into 
it (IMoriva et al.ll2010t ITanaka et aLl l2012c iMo riva et all 
I2013D (see also iTominaga et al.ll201ll) . 



We have modeled Pop III Type Iln SNe and their 
light curves and spectra with the Los Alamos RAGE and 
SPECTRUM codes in order to calculate their detection 
limits in rcdshift and their NIR signatures. In Section 
2 we describe our explosion and shell models and how 
they are evolved in RAGE. Blast profiles, light curves 
and spectra are examined in Section 3, and in Section 4 
we compare light curves from our models with those of 
Type Iln SN candidates observed in the local universe. 
In Section 5 we calculate NIR light curves for Pop III 
Type line at high z and determine their detection limits 
as a function of redshift. In Section 6 we conclude. 

2. NUMERICAL MODELS 

We take the z40G SN from IWhalen et all (|2012cf ) 
(hereafter WET12) to be our fiducial explosion model. 
It is a zero-metallicity, 40 Mq, 2.4 x 10 51 erg CC SN. 
This progenitor was evolved from the beginning of the 
main sequ ence up to the poin t of explosion in the Ke- 
pler code ([Weaver et al.l 11978c IWooslev et all 120021 ). at 
which point the SN was artificially triggered and fol- 
lowed until the end of all nuclear burning. Profiles for 
the blast, which at this point was still deep in the star, 
were then mapped onto a 2D grid i n the CASTRO adap- 
tive mesh refinement (AMR) code (jAlmgren et alll2010D 
and evolved until just before shock breakout to capture 
mixing inside the star due to Rayleigh- Taylor instabili- 
ties. We begin this study by spherically averaging the 
final CASTRO density, energy, velocity and mass frac- 
tions and mapping them onto a ID spherical AMR grid 
in RAGE together with the surrounding star, its wind 
and a variety of dense shells. 

We adopted z40G, a red supergiant progenitor rather 
than a blue compact giant star like those thought to be 
the origin of LBV outbursts, because its larger surface 
area gives an upper limit to the luminosity from shock 
breakout and hence the number of photons that initially 
escape through the shell. The dynamics of the collision 
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with the shell is not sensitive to the structure of the star. 
Had we instead used the u40G explosion a greater frac- 
tion of the breakout transient would have gotten through 
the shell because more of it would have been x-rays 
(shocks breaking out of compact stars are hotter than 
those breaking free from red giants for a given explosion 
energy). On the other hand, blue stars have smaller sur- 
face areas and this compensates for higher shock temper- 
atures in the flux. None of these issues impact detection 
limits because the breakout pulse is completely absorbed 
by the neutral IGM at high redshift. 

2.1. RAGE 

RAGE (Radiati on Adaptive Grid Eulerian; 
iGittings et all [20081 ) is a multidimensional adaptive 
mesh refinement radiation hydrodynamics code de- 
veloped at Los Alamos National Laboratory (LANL). 
It couples second order conservative Godunov hydro- 
dynamics to grey or multigroup flux-limited diffusion 
(FLD) to model strongly radiating flows. Our RAGE 
root grid has 200,000 zones with a resolution of 4.0 x 
10 10 cm and reflecting and outflow conditions on the 
inner and outer boundaries, respectively. Our choice of 
mesh ensures that all important features of the shock, 
the star, and the shell are resolved by at least 10 zones 
and that 5000 zones are allocated from the center of 
the grid to the outer edge of the shock at setup. As in 
WET12, we allow up to five levels of refinement so each 
feature can be further resolved by up to 320 additional 
zones if necessary, and we again periodically resamplc 
the explosion onto larger grids to accommodate its ex- 
pansion and speed up the calculation. All physics in the 
WET12 models are used here: multi-species advection, 
grey flux- limited diffusion radiation transport with Lo s 
Alamos OPLIB atomic opacitie£E (jMagee et all 119951 ). 
2-temperature physics, and energy deposition due to 
radioactive decay of 56 Ni. 

We evolve the explosion through breakout from the 
star, collision with and propagation through the shell, 
and expansion into the IGM out to 500 days. When we 
post-process dumps from RAGE with SPECTRUM to 
obtain light curves, we sample shock breakout from the 
star with 50 spectra at evenly spaced times that bracket 
the thermal transient and 200 spectra at logarithmically 
spaced times out to 500 days. Computing spectra in 
this manner allows us to model with detailed opacities 
how the shell attenuates both the initial breakout pulse 
and the flash from its own lower layers when the ejecta 
collides with it. 

2.2. SPECTRUM 

RAGE pro files are importe d into the LANL SPEC- 
TRUM code (jFrev et al.ll2013l ) in order to calculate spec- 
tra with monochromatic OPLIB atomic opacities that 
capture detailed emission a n d abs orption line structure. 
As explained in iFrev et all (|2013l ). densities, velocities, 
and gas and radiation energy densities from the most 
refined levels in the RAGE AMR hierarchy are first ex- 
tracted and reordered by radius into separate data files. 
These profiles are then mapped onto a 2D grid in radius 
and /j, = cos 9 in SPECTRUM. The spectrum calculation 

11 http:/ /aphysics2/www. t4.lanl.gov/cgi-bin/opacity/tops. pi 




1 (T 5 L I 

10 13 10 14 10 15 10 16 10 17 
r (cm) 



Fig. 1. — Density profiles for the shells considered in our study 
(models A00 through A04). The SN shock is visible at 4.0el3 cm 
and the surface the 40 Mq star is visible at 1.3el4 cm. 



TABLE 1 

POP III ClRCUMSTELLAR SHELL MODELS 



model 


Msh {Mq) 


rh w (Mq 1 yr) 


v w (km s x ) 


A00 


0.1 


10" 4 


200 


A01 


1 


10~ 4 


200 


A02 


6 


10~ 4 


200 


A03 


10 


10~ 4 


200 


A04 


20 


10~ 4 


200 



is done with 5200 radial bins and 160 angular bins, the 
same number as in WET12 but with a slightly modified 
gridding strategy. The simulation volume is divided into 
3 regions: inside the r = 20 surface, between the r = 20 
surface and the radiation front, and beyond the radiation 
front. The position of the radiation front is taken to be 
the outermost cell with a temperature above 0.0292 eV, 
approximately three times the background temperature 
of 0.01 eV (136 K). The radius of the r = 20 surface is 
calculated from the outer edge of the grid assuming k = 
0.2 cm 2 /g, that due to electron scattering in gas at pri- 
mordial composition, 76% H and 24% He by mass. To 
avoid allocating to many zones to the cold dense shell 
and under- resolving the shock prior to its impact with 
the shell, the r = 20 surface calculation excludes the 
density of the shell before the collision. If a zone in the 
shell has not received additional momentum from the 
shock then the density in that cell is replaced with what 
the wind density would be at that location. This allows 
SPECTRUM to adequately resolve the photosphere of 
the ejecta prior to its collision with the shell and then 
the shell itself when it becomes the dominant source of 
flux. 

2.3. Shell Structure 

Pop III stars are not normally thought to lose much 
mass over their lifetimes because there a re no line driven 
winds in their metal- free atmospheres (jKudritzkil 120001 : 



Vink et al.1 [200lt iB araffc c t~aTl [200lt iKrticka k Kubatl 
20061: lEkstrom et all I2008H . However, these studies do 
not exclude the possibility of violent pulsational mass loss 
late in the life of massive primordial st ars. For simplic- 
ity, we adopt the shell structure used in Ivan Marie et aTl 
(|2010f K a transient, high- mass flux that interrupts the 
more diffuse wind blown by the star before and after the 
ejection. The two winds have the same constant velocity 
but different uniform mass loss rates that together cre- 
ate a density profile that is a simple superposition of two 
winds: 



p(r) 



. m j" — if r < n and r > 



m ^ if ri < r < r 2 . 



4irr' 2 v u 



Here, m w and rh s h are the mass loss rates of the wind 
and the shell, where 



M. 



sir 



dt. 



(1) 



and v w is the wind speed. The two radii r\ and r 2 mark 
the inner and outer surfaces of the shell, 



^1 V w t en d 

r2=v w (t end + dt sh ), 



(2) 
(3) 



where t en d is the time between the end of the ejection 
and t he SN, and dt s h is th e duration of the ejection. 
As in Ivan Marie et all (|2010f >. t end and dt sh are 2 yr in 
all our models and we vary the density of the shells by 
adjusting M s h- We take the shells to be primordial, 76% 
Hand 26% He by mass fraction. Density profiles for the 
5 shells in our study together with the 40 Mq star are 
shown in Figure [TJ and we summarize the properties of 
the shells in Table [Q 

The aim of our numerical campaign is to explore the 
observational signatures and detection thresholds of Pop 
III Type Iln SNe, not perform an exhaustive survey of 
such explosions. We adopt t his type of shell to com pare 
our light curves with those of lvan Marie et al.l (|2010h and 
because it is similar in mass and structure to the shell in- 
ferred to exist around n Carinae from observations. Den- 
sity profiles for actual LBV eruptions are likely more 
complicated, with fast winds preceding and following 
much slower outbursts that create multipl e shocks and 
rarefa ction zones similar to those found in iMesler et al.1 
(|2012f) . Radiative cooling would also flatten the shell into 
a colder and denser structure than the ones shown here 
if dust and metals are present. Type Iln SNe in the local 
universe will be examined in a forthcoming study. 

2.4. Ionization State of the Shell 

Before launching our runs in RAGE we performed a 
separate test to determin e if UV radiation fro m the star 
ionizes the shell (see, e.g.. lWhalen et aLl l2004f). The ion- 
ization state of the shell determines its opacity to the 
SN before and after its collision with the shell. It may 
also influence how efficiently the kinetic energy of the 
ejecta is transformed into radiation upon impact with 
the shell, and hence its luminosity. Both issues are rele- 
vant because the progenitor can be tens of solar masses 



and therefore a very luminous source of ionizing UV pho- 
tons, and there is little if any dust in the vicinity of the 
star to attenuate them. 

We use the ZEUS-MP code to determine if th e star ion- 
izes its shell (jWhalen k Normadl2lM l2008al lbT). ZEUS- 
MP self-consistcntly solves hydrodynamics, noncquilib- 
rium H and He chemistry and ionizing UV transport to 
propagate cosmological ionization fronts. We consider 
a 250 M Q Pop III star in shell A00 from Table ffl The 
wind and shell are initialized on a ID spherical mesh with 
200 uniform zones and inner and outer boundaries at 1.3 
x 10 13 cm (the surface of the star) and 3.0 x 10 15 cm 
(the outer surface of the shell 2 yr after the end of the 
ejection). We use multifrequency UV transport, with 40 
bins uniformly partitioned in energy from 0.255 to 13.6 
eV and 80 bins logarithmically spaced from 13.6 to 90 
cV. The blackbody spectrum of the star is normalized 
to ionizing photon em ission rat e s, sur face temperatures, 
and luminosities from iSchaererl (|2002f ) . The shell is illu- 
minated for 4 yr, the time from the onset of ejection to 
the SN. This treatment is approximate, given that the 
shell is initially closer to the star and exposed to higher 
fluxes just after expulsion. Because a 250 Mq star is 
more luminous than a Type Iln progenitor, its ionizing 
flux is the extreme upper limit that could be reasonably 
applied to the shell. 

We find that the radiation front from the star easily 
ionizes the wind on timescales of a few hours but is halted 
by the shell without ionizing even one zone of it. Since 
the star is incapable of ionizing the least massive shell in 
our study, we take all the shells in our RAGE models to 
be neutral. We note that had the shells been ionized they 
would fully recombine by the time the explosion reaches 
them. This is evident from the recombination timescales 
in the gas, 

1 



n e a(T) ' 



(4) 



where 



a(T) = 2.59 x W'^T 



-13^-0.75 cm -3 s -l 



(5) 

is the case B recombination rate coefficient for hydrogen, 
n e is the electron number density and T4 is the temper- 
ature in units of 10 4 K. With densities of 10 9 cm~ 3 in 
the A00 shell and ionized gas temperatures T ~ 10,000 
K, t rec is about an hour. 

3. POP III UN SN LIGHT CURVES AND SPECTRA 

We show bolomctric light curves for all five explosions 
in the source frame out to 60 days and 500 days in the 
left and right panels of Figure [21 respectively, and spec- 
tra for shock breakout from the surface of the star in 
Figure [3l Blast profiles for the A04 explosion are shown 
in Figures [4] - [6] The evolution of the SN can be par- 
titioned into three distinct phases: (1) the collision of 
a radiative precursor with the shell; (2) the collision of 
the ejecta with the shell and its propagation through the 
shell; and (3), breakout from the outer surface of the 
shell and expansion into the surrounding medium. 

3.1. The Radiative Precursor 

As shown in the left panel of Figure [51 some flux 
from shock breakout from the surface of the star passes 
through all five shells. The breakout transient, shown at 
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Fig. 2. — Bolometric luminosities for all 5 Type Iln explosions together with the z40G explosion. Left panel: out to 60 days. Right panel: 
out to 500 days. Shock breakout from the star itself appears as the luminosity spike at 4.5 days, the collision of the radiative precursor 
with the inner surface of the shell is visible as the small uptick in luminosity at 9 days, and the collision of the cjccta with the shell causes 
the jump in luminosity at 30 - 40 days. Breakout of ejecta from the outer surface of the shell is visible as the surge in bolometric luminosity 
at 60 - 125 days on the right. 

precursor, the feature at 4.0 x 10 14 cm at 5.6 days in the 
density and velocity profiles in Figure [4l reaches with 
the inner surface of the shell at 9 days, well before the 
ejecta. The collision creates a temperature spike at the 
inner surface of the shell that is manifest as the small 
bump in luminosity just after the breakout transient in 
all five shell light curves. After the initial collision, the 
wispy outermost layers of the star continue to pile up in 
a thin layer at the inner surface of the shell, forming a 
reverse shock that reverberates back and forth through 
the layer. This can be seen in the temperature profile 
across the layer: the small spike in temperature at 1.3 
x 10 15 cm, the upper surface of the layer, at 9.0 days 
rebounds to 9.0 x 10 14 cm, the lower surface, by 10.9 
days. Slightly fewer high-energy photons get through 
the shell at 10.9 days than 9 days because the reverse 
shock is at the bottom of the thin layer and more of its 
radiation is absorbed. 

The reverberation of the shock back and forth in this 
layer causes the flickering in the light curve from 8-32 
days. The luminosities arc relatively steady over these 
times because the thin layer is trapped at the inner sur- 
face of the shell. However, luminosities for less massive 
shells are higher because they allow more radiation to 
pass. The ripples in the A02 - A04 light curves are ab- 
sent from the A00 and A01 light curves for two reasons. 
First, as we discuss in greater detail below, the ejecta 
never overtakes the inner surface of the A00 shell because 
its radiation displaces the low-mass shell outward. Sec- 
ond, the radiation front propagates to greater distances 
through the shell because of its lower density, so there 
is a larger separation between the front and the r = 20 
surface. As a result, the SPECTRUM code may not re- 
solve the radiating region of the flow and capture small- 
amplitude fluctuations in luminosity in diffuse shells as 
well as in more massive shells. In the A02 - A04 shells the 
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Fig. 3. — Spectra at shock breakout from the surface of the star 
(4.5 days) for all five Type Iln explosions. The absorption features 
at 1000 - 3000 A due to the shell cause the drop in bolometric 
luminosity at 4.5 days. 

5 days in the left panel of Figure [21 is ~ 5 x 10 45 erg 
s _1 but is reduced to 10 43 erg s _1 by the A00 shell and 
11 erg s _1 by the A04 shell. This attenuated flux is 



to 10 4 

unlikely to be detected because it is relatively dim and 
because the shell becomes much brighter later on. As we 
show in Figure [21 most of the attenuation is due to the 
absorption of 1000 - 3000 A photons. The bolometric lu- 
minosity of the shock falls to ~ 10 42 ' 5 erg s _1 after about 
a day and decays slowly thereafter. 

The photon pulse blows off the outer layers of the star 
after breaking free of the shock and drives them across 
the gap between the star and the shell. This radiative 
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Fig. 4. — The collision of the radiative precursor, the outermost layers of the star blown off by radiation at shock breakout, with the 
inner surface of the A04 shell. Upper left panel: spectra; upper right panel: densities; lower left panel: radiation temperatures; lower right 
panel: gas velocities. As the precursor piles up against the inner surface of the shell, a shock (visible as the peaks in radiation temperature 
at 9.0el4 cm at 9.0 days and 1.3el5 cm at 10.9 days) reverberates back and forth through the accumulated gas. 



radiation front and the r = 20 surface are much closer 
together and the radiating region is better resolved. 

3.2. Collision with the Shell 

We show the ejecta just before and after its collision 
with the shell at 27 days and 49.2 days in Figure [5] When 
it reaches the inner surface, the ejecta drives a strong 
shock into the shell, heating it to ~ 7.5 eV. The shock 
drives the inner surface outward, piling it up into the 
density spike visible at 1.5 x 10 15 cm at 49.2 days. This 
7.5 eV spike is the source of the jump in luminosity at 
32 days in the left panel of Figure [5J The magnitude of 
the jump depends on the density of the shell, with diffuse 
shells allowing more radiation to pass through them. The 



light curves appear to plateau from 32 - 60 days, but the 
luminosity rapidly rises as the shock propagates through 
the shell and there is less gas between the shock and the 
outer surface of the shell to absorb or scatter the photons, 
as we show at 32 - 100 days in the right panel of Figure^ 
The rise is most pronounced in the densest shells because 
they initially attenuate the most radiation. 

As the ejecta plows up the shell a reverse shock forms 
and detaches from the forward shock, backstepping into 
the interior in the frame of the flow. The reverse shock 
can be seen at 2.0 x 10 15 cm in the density and tempera- 
ture profiles at 93.5 days. As it recedes from the forward 
shock, the reverse shock loses pressure support because 
of radiative cooling by H and He lines in the shocked gas 
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Fig. 5. — The collision of the ejecta with the A04 shell. Upper left panel: spectra; upper right panel: densities; lower left panel: radiation 
temperatures; lower right panel: gas velocities. The formation of a reverse shock as the ejecta plows through the shell can be seen in the 
double peak in the velocity profiles at 49.2 and 93.5 days. 



and it retreats back toward the forward shock. However, 
as the forward shock continues to plow up the shell the 
reverse shock is again driven back into the interior. The 
cyclical heating and cooling of postshock gas associated 
with the oscillation of the reverse shock causes the rip- 
ples in the luminosities from 32 - 100 days in the A02 - 
A04 light curves. The pe riod of oscillation is g overned by 
cooling rates in the gas l|Imamura et al.lll984r ) and is in- 
dependent of shell mass, but the amplitude of oscillation 
is somewh at correlated with sh ell density (see also sec- 
tion 4.1 of lAnninos et al.lll997l ). Such features are also 
found in Lyman alpha emission by primordial SN rem- 
nants at later times as they sweep up neutral gas in halo s 
on larger scales (see Figure 11 in iWhalen et al.ll2008d ). 
The stability of radiative shocks in astrophysical contexts 



has been well studied (e.g.. iChevalier fc Imamuralll982[ ). 
We note that the spectrum becomes harder from 27 to 
93.5 days even though the shock cools because there is 
less absorption by the shell as the shock plows through 
it. 

3.3. Breakout from the Shell 

The shock breaks free from the outer surface of the 
shell at times that depend on the mass of the shell, from 
~ 45 days in A01 to ~ 120 days in A04. As shown in 
the velocity profiles at 120 and 139 days in FigureEJ the 
shock abruptly accelerates in the sharp density drop just 
outside the shell. Photons from the shock also stream 
into the surrounding medium, creating the second jump 
in luminosity to ~ 10 43 5 erg s _1 that persists for 50 - 100 
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Fig. 6. — Shock breakout from the A04 shell. Upper left panel: spectra; upper right panel: densities; lower left panel: radiation 
temperatures; lower right panel: gas velocities. The initial trapping and subsequent escape of radiation into the IGM can be seen in the 
absence of a plateau in radiation temperature ahead of the shock at 129 days and its appearance and propagation at 139 and 154 days. 
Breakout from the shell is most evident in the jump in gas velocities between 120 and 139 days. 

154 to 500 days both of these effects are evident in the 
spectrum. The peak of the spectrum shifts to longer 
wavelengths but more and more high energy photons es- 
cape into the IGM. At 500 days the spectrum is nearly 
blackbody when its cutoff at high energies has fallen be- 
low the Lyman limit and ambient densities are very low. 
From 250 - 500 days all five light curves gradually dim 
as the shock expands and cools. 



days in the right panel of Figure [5] Because the peaks 
coincide with shock breakout, they also occur at 45 - 
120 days. Diffuse shells result in broader peaks because 
they allow photons to escape from greater depths and 
earlier times in the shell. As we show at 120 days in the 
temperature profile, low-energy photons begin to leak 
through the shell well before shock breakout at 139 days. 

Radiation from the shock again blows off the outer lay- 
ers of the shell, creating the density peak at the edge of 
the radiation front at 4.0 x 10 15 cm at 139 days. These 
photons ionize the envelope beyond the shell, allowing 
more energetic photons with A > 100 A to pass through, 
as shown in the spectrum at 154 days. As the shock 
expands it cools, and the ambient density falls. From 



3.4. A00 

The A00 explosion evolves somewhat differently be- 
cause its shell is so diffuse. The radiative precursor 
crosses the gap between the ejecta and the shell in ~ 
9 days, but radiation from the ejecta has already perme- 
ated and warmed the shell. The precursor piles up in a 
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Fig. 7. — The collision of the ejecta with the A00 shell. Upper left panel: spectra; upper right panel: densities; lower left panel: radiation 
temperatures; lower right panel: gas velocities. Note from the profiles that the ejecta does not fully close the gap with the shell; radiation 
in the gap instead propels the shell forward. 



thin layer at the inner surface of the shell as in the other 
cases, but as the ejecta approaches this surface the gap 
is never fully closed. Radiation in the gap both displaces 
the inner surface outward and drives a reverse shock into 
the ejecta, as we show at 20.3, 40.9 and 59.3 days in 
Figure [7j The ejecta drives the inner surface forward 
while remaining separated from it by radiation forces. 
As shown at 59.3 days, there is still a gap between the 
ejecta and the shell, and radiation from the shock has 
leaked through the shell. By this time the outer edge of 
the shell has also begun to be displaced outward and it 
is quickly accelerated to the same velocity as the ejecta. 

This sequence of flow has several consequences for the 
A00 light curve. First, as seen in FigureEJ the luminosity 
from 8-32 days is higher because more radiation from 



the collision of the precursor with the shell gets through 
the shell. Next, because the ejecta following the precur- 
sor is prevented from colliding with the shell by radiation 
in the gap between them, the rise in luminosity as the 
ejecta approaches the shell is very gradual, without the 
well-defined jump at ~ 32 days in the A02 - A04 light 
curves. Besides being more transparent to radiation, the 
A00 shell is also swept up more quickly than the others. 
The shock that is driven by the radiative gap reaches the 
outer surface of the shell at 60 days, when its luminosity 
peaks. From 32 - 60 days the shock is far less lumi- 
nous than in more massive shells because the A00 shell 
is being rapidly accelerated to the velocity of the ejecta. 
Consequently, the kinetic energy of the ejecta is not ef- 
ficiently converted into thermal energy, and the shock is 
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weaker and dimmer while it is in the shell. Finally, the 
luminosity jump at breakout in the other shells is absent 
in the A00 light curve because the ejecta never actually 
passes through the shell. Supcrnovae in low-mass shells 
are much dimmer than those in shells above 1 Mq. 

The A00 light curve rcbrightcns between 110 and 220 
days because the photosphere sinks into the ejecta and 
encounters layers of higher density and temperature, so 
the broad peak at 220 days is due purely to optical depth. 
A similar feature is visible at 240 days in the A01 light 
curve. The evolution of the A01 explosion is intermediate 
to that of A00 and the others. Its luminosity gradually 
rises from 8-32 days because radiation from the impact 
of the outer layers of the star with the shell gets through 
the shell and intensifies as gas piles up at the inner sur- 
face. Next, because a radiativcly supported gap again 
softens the collision of the ejecta with the shell, the rise 
in luminosity upon impact is again more gradual than 
with more massive shells but more prominent than in 
the A00 light curve. For shell masses greater than 1 Mq , 
breakout luminosities from the shell are 10 43 5 - 10 43 ' 8 
erg s _1 . 

3.5. Earlier Simulations 

Not surprisingly, radiation transport leads to impor- 
tant qu alitative differenc e s betw een our light curves and 
those of Ivan Marie et all (120101) . even though we use the 
same shells in our models. Tvan Marie et al.l(|2010D equate 
energy losses due to optically thin radiative cooling in the 
gas with the bolometric luminosity of the SN remnant 
and hence do not account for absorption by the optically 
thick shell. These losses are greatest when the gas is 
most strongly shocked, when the ejecta strikes the inner 
surfac e of the shell. Consequently, the Ivan Marie et al.l 
pOlOh light curves peak upon initial collision when in 
reality most of these photons are absorbed in the lower 
layers of the shell (compare their Figure 11 at 25 days 
to the right panel of Figure [2] at 32 days). On the 
ot her hand, because colli sional gas cooling scales as p 2 
in Ivan Marie et al.1 (|2010T l , their bolometric luminosities 
plummet when the shock breaks out into the low-density 
medium beyond the shell (compare Figure 11 at 60 - 
100 days to Figure [5] at 50 - 100 days). As shown in 
Figure [21 the light curve actually peaks at this point be- 
cause photons previously trapped in the shock are free 
to stream into the IGM, and as the shock rushes down 
the density gradient it heats and radiates additional en- 
ergy. When the SN collides with the shell our light 
curves are therefore initially dim and later peak upon 
breakout from the shel l, while the converse is true of the 
van Marie et al.l (|2010l) light curves. We note that our 
shock prop agates through the she lls in about the same 
times as in Ivan Marie et aTl (|2010D . 

Interestingly, both approaches yield similar widths for 
the light curves at peak luminosity. The width of the 
peak is closely tied to the prop agation time through the 
shell in Ivan Marie et aTl (|2010f ) but is largely a function 
of the subsequent expansion and cooling of the remnant 
in our models. Since shock temperatures peak at only 
~ 8 eV upon impact with the shell, we do not find any 
x-ray emission after breakout from the surface of the star 
but the SN is bluer out to much later times than in z40G 
by itself. This property makes Type Iln SNe ideal candi- 
dates for high redshift detection. We note that some of 



the lig ht curves in Figures 10 and 11 of Ivan Marie et aTl 
((2010ft exhibit the same flickering as in ours, probably 
also because of oscillations of reverse shocks due to ra- 
diative cooling. They are especially prominent as the 
shock plows through the A00 shell in Figure 11, when a 
reverse shock is l ikely to detach from the forward shock. 
Finally, because Ivan Marie et al.l ((2010ft model the SN 
profile as a free expansion and do not consider breakout 
from the star, none of the features from 8-32 days in 
our light curves are present in theirs. 

We note that while the shock never reaches tempera- 
tures at which it would emit x-rays in our models, this 
in part is due to the fact that we assume that ions and 
electrons are closely coupled. Because the electrons and 
the photons are also closely coupled, the ions rapidly lose 
energy to the radiation field. Because e gas = CyT and 
e-rad = aT 4 , large transfers of energy into the radiation 
field do not result in large changes in its temperature, 
so the radiation and matter temperatures remain mod- 
est and there are only mild deviations between them in 
our simulations. In reality, the ions may not couple effi- 
ciently to electrons when the shock collides with the inner 
surface of the shell or breaks free from its outer surface 
and the shock may heat them to much higher tempera- 
tures than the radiation temperature. Under these cir- 
cumstances the ions may drive a small fraction of the 
electrons out of equilibrium with the radiation field and 
cause them to emit x-rays in real SN remnants. 

Comparison of the more recent light curves of 
iMoriva et aTl (12013ft with ours is more problematic, in 
part because they assume much higher explosion ener- 
gies (10 - 50 x 10 51 erg, rivaling those of pair-instability 
SNe) and because they adopt a different structure for 
the shell. As a result, their peak luminosities are 5 - 
10 times higher than ours but have similar widths. Our 
bolometric light curves are otherwise qualitatively simi- 
lar to theirs (compare their Figure 5 to the left panel of 
Figure [2]) . Their large explosion energies are necessary 
to achieve the peak bolom etric luminosities of SN 2006gy 
(jSmith et al.ll2007l l2008bft and they are not normal core- 
collapse events. 

4. TYPE UN SNE IN THE LOCAL UNIVERSE 

We now compare bolometric luminosities for our 
five explosions to those of Type I ln SNe discovered 
in the local univers e: SN2006gy (ISmith et al.1 120071 
I2008bft. SN2006tf (ISmith et all I2008al). SN2007pk 
dPritchard et al.l |2012|), SN2008am dChatzopoulos et al.l 
12011ft. SN2 10al (jCooke et al.l 12010ft. SN2010jp 
(ISmith et all 120121). and SN2011ht (jRoming et al.l 
120121 Humphreys et alJl2012t) . SN2006gy and SN2006tf 
are the most luminous explosions in this class while the 
others exhibit more typical peak luminosities. We group 
luminosities for A00 and A01 with those of normal Type 
Iln S Ne in the left panel of Figure |5] (see iKiewe et al.l 
120121 for additional examples of Type Iln SN light 
curves) and light curves for A02, A03 and A04 with 
those for the superluminous SN2006gy and SN2006tf 
in the right panel of Figure |8] The rise of the light 
curve is evident in SN2006gy and SN2011ht; the other 
five datasets only show its decline, making it difficult 
to pinpoint their explosion times. We therefore assign 
them times that best align them with our simulations. 
The observations of the normal brightness Type Iln SNe 
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Fig. 8. — Bolomctric light curves for Type Iln SN candidates in the local universe overlaid with the A00 - A04 light curves. Left: 
less luminous Type line and the A00 and A01 models. Right: superluminous Type Iln candidates with the A02 - A04 models. In the 
observations, there is a hint of shock breakout from the surface of the star with SN2011ht at 10 days that is absent in the other data. 
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Fig. 9. — Observations of z ~ 2.2 Type Iln SN candidates vs. simulations. Left: r band. Right: i band. In both panels, the first peak is 
due to the collision of the ejecta with the inner surface of the shell and the second peak is due to shock breakout from the shell. The first 
peak is mostly absent in the r band because the photons are blueward of the 3000 A absorption feature of the spectra in Figure [3] and are 
cut off by the shell. 



either have very few data points (SN2008am) or cover 
only short periods of time (SN2007pk, SN2010al, and 
SN2010jp). 

The central engines of Type Iln SNe, the core- 
collapse explosion, do not vary much with met allicity 
(IChieffi fc Lin7o7gil[200l IWooslev fc Hegerl I2007I) (note 
also Figure 1 of IWhalen fc Frverll2012f ). It is primarily 
the structure of the shell and its cooling properties that 
differentiate Pop III Type Iln supernovae from those to- 
day. Fine-structure cooling by metals flattens shells into 
colder, thinner and denser structures than do H and He 
lines. The thickness of the shell governs the width of the 
main luminosity peak, so Type Iln events today would 
likely exhibit narrower peaks and sharper declines in 
bolometric luminosity, as with SN2010jp and SN2010al. 



On the other hand, inefficient gas cooling keeps the shock 
hot and bright, with much slower decays in luminosity at 
later times like those in Figure [8] More realistic treat- 
ments of the ejection would impose additional features 
on both zero-metallicity and enriched shells that are not 
present in our models. Slow outbursts are usually pre- 
ceded and followed by much faster winds. The wind in 
front of the shell detaches from and races ahead of it, 
creating a rarefaction zone, while the wind behind the 
shell piles up at it s inner surface and forms a hot termi- 
nation shock (see IMesler et al.l |2012j ) . These structures 
will imprint additional features on the light curves that 
have not been captured by any simulations to date. Type 
Iln SNe in more realistic shells in the local universe will 
be pursued in future simulations. 
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In spite of these limitations, the A00 and A01 light 
curves are consistent with those of the five less luminous 
Type line. They fall between A00 and A01, with rates 
of decline that are similar to those in the simulations. In 
particular, A01 is an excellent match to SN2007pk. Be- 
cause SN2011ht was detected during its rise and observed 
for longer times, it places tighter constraints on our mod- 
els. Its luminosity peak falls almost directly between 
those of A00 and A01 on a log scale, and it has about 
the same width. However, the ratio of peak to plateau 
luminosities is smaller for SN2011ht than for the two sim- 
ulations. The z40G explosion in a shell with a mass of ~ 
0.4 M Q would yield the best match to SN2011ht. 

SN2006gy, another Type Iln SN that was observed dur- 
ing its rise to peak luminosity, is nearly a factor of 3 
brighter than A04, our brightest light curve. As noted 
earlier, bolometric luminosities rose only slightly as the 
mass of the shell went from 6 to 20 Mq, so more mas- 
sive shells will not yield better agreement with SN2006gy. 
The dip at the beginning of the SN2006gy light curve 
may be photons from shock breakout from the surface of 
the star filtering through the shell, which only happens 
with less massive shells. Taken together, these two facts 
suggest that a more powerful SN is needed to explain 
SN2 006gy, not a more m assive shell (which is consistent 
with lMoriva et al.ll2013l ). SN2006tf, the other superlumi- 
nous explosion, is only marginally brighter than A04. Its 
bolometric luminosities are close to those of A04 at early 
and late times, but A04 falls more rapidly and then en- 
ters a plateau while SN2006tf declines more steadily. As 
mentioned above, the plateau is likely due to inefficient 
H and He cooling in the A04 shell and would probably 
disappear if t he gas was enrich e d wit h metals. The A04 
light curve in Ivan Marie et all (pOlOf ) is in basic agree- 
ment with SN2006gy for shell masses of 20-24 M and 
wind velocities of 190 km s _1 . 

Our light curves are in general agreement with recent 
Type Iln SNe, and some observations match our simu- 
lations extremely well. Nevertheless, we do not expect 
exact agreement because we only used one explosion and 
shell, and varied only the density of the shell. The shells 
in our models are also somewhat different from those in 
real explosions. From Figure [5] it is clear that many of 
the properties of the shell can be extracted from the light 
curves, so they can be powerful probes of the circumstel- 
lar environment of the explosion. 

5. Z ~ 2 TYPE UN SNE 

In Figure [5] we compare r-band and i-band light 
curves for the A01 - A04 runs with those of Type Iln 
SNe recently discovered at 1.9 < z < 2.4: SN 0224- 
0457, SN 0224-0426, SN 1420+5252, and SN 2214-1807 
(jCooke et alJl2009L 12012ft with the Low Resolution Imag- 
ing Spectrometer (lOke et all 119951: iSteidel et all 120041 : 



McCarthy etHI 119981 IRockosi et all 12010ft . These SNe 
reach peak AB magnitudes of 24.5 - 26 and are visible 
for 200 - 250 days in both bands. Our simulated r-band 
light curves exhibit two peaks: a brief initial peak lasting 
no more than 50 days with magnitudes below 28, and a 
second brighter and longer peak that reaches magnitudes 
of 27 - 28 and lasts up to 500 days. The i-band light 
curves also exhibit two peaks but the first one is longer 
and brighter than in the r-band: 50 - 150 days with mag- 
nitudes of 26.5 to 27.5. In both bands the first peak is 



due to the collision of the SN ejecta with the inner edge 
of the shell and the second peak is due to shock break- 
out from the shell, when its photons are suddenly able to 
stream freely in the low density wind. The first peak in 
the i band becomes brighter and longer as the shell mass 
decreases. In Figure [5j the shock reaches temperatures 
of ~ 10 eV and becomes quite luminous upon its collision 
with and initial advance into the shell. However, these 
photons must filter through the shell, and the number 
that escape depends strongly on wavelength. As shown 
in Figure[51 the opacity of the shell imposes a sharp cutoff 
on the spectrum at ~ 3000 A. This explains the promi- 
nence of the first peak in the i band. The photons that 
are redshifted into the i band from z — 2.2 originate from 
the brightest region of the spectrum just redward of the 
3000 A cutoff. There is little luminosity in the r band 
because they are blueward of the cutoff and absorbed by 
the shell. 

The first peak in the synthetic i-band light curves bet- 
ter matches the observations than the longer and dimmer 
late time peak that appears in both bands. If these tran- 
sients arc Type Iln events it is likely that the initial colli- 
sion with the shell is being observed, and breakout from 
the shell would be seen later if the surveys were extended. 
We note that at z = 2.2 the progenitors would be Pop 
II stars, not pristine Pop III stars, and the presence of 
even small amounts of dust or metals in the shell might 
have large effects on the opacity and further reduce the 
magnitude of the first peak in the i-band. In addition, 
our simulation suite only explored the effect of changing 
the mass in simple analytic shells. The light curves also 
depend on the energy of the SN explosion, the distance 
from the SN to the shell, and the thickness of the shell. 
Therefore, it is not surprising that our 5 models do not 
exactly match the observational data. Future simula- 
tions are planned to explore a larger parameter space for 
shell collision SNe in the local universe. 

6. NIR LIGHT CURVES 

We calculate NIR light curves for Pop III Type Iln 
SNe with the synthetic photometry code described in 
ISu et all (|2011[ ). We redshift each spectrum to the de- 
sired z before removing the flux absorbed by i nterven- 
ing ne utral hydrogen clouds with the method of iMadaul 
(|1995|) . Each spectrum is then dimmed by the required 
cosmological factors. Our algorithm linearly interpolates 
the least sampled data between the input spectrum and 
filter curve. It has additional capabilities such as redden- 
ing by dust that are not used here. 

We show NIR light curves for all 5 explosions at z = 7, 
10, 15, 20 and 30 in Figures [T0l and [TT1 In each case we 
plot the NIR signal in the optimum filter for its detection, 
which in all cases is redward of the Lyman limit at that 
redshift. At z = 7 and 10, two peaks are visible in the 
A01 - A04 explosions, a narrow, brighter peak at 500 - 
1000 days and a dimmer, broader peak at 2000 - 2500 
days. The first is due to shock breakout from the outer 
surface of the shell. The second occurs as the redshifted 
spectral peak evolves downward in wavelength through 
the NIR as the shell later expands and cools. The first 
peak is present but not visible at higher redshifts. Its 
luminosity does not change with shell mass, but more 
diffuse shells have broader peaks because photons from 
the shock can escape such shells before the shock, with 
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Fig. 10. — NIR light curves for all five shell explosions. Upper left panel: z = 7; upper right panel: z = 10; lower left panel: z = 15; 
lower right panel: z = 20. The dashed line at AB mag 32 indicates the JWST NIRCam photometry limit. The first peak at z = 7 and 
10 is due to the initial collision with the shell and cannot be seen at z > 10. The second peak is caused by breakout from the shell and is 
visible out to z > 20. 



lead times that are inversely proportional to the density 
of the shell. The first peak appears at earlier times with 
low-mass shells because the shock breaks free of them 
sooner. The A00 light curve does not exhibit this peak 
because there is no breakout from the shell. In contrast, 
the luminosity of the second peak rises with shell mass 
but its width is relatively uniform. As expected, this 
second peak occurs at later times at higher redshifts. 

With NIRCam photometry limits of AB magnitude 31 
- 32, JWST will be able to detect all five explosions out 
to z ~ 20. With proposed NIR detection limits of AB 
magnitude 26.5 - 27, WFIRST and WISH will only be 
able to observe such events out to z ~ 7 - 10. At z = 7 
and 10 the first peak is visible for 250 - 500 days and the 



second peak is visible for ~ 1500 days. At z = 15 and 20 
the second peak can be seen for ~ 2000 days but is about 
a magnitude dimmer. The peaks rise as quickly as they 
fall, with durations that are comparable to likely proto- 
galactic survey times of 1 - 5 yr. Because the spectrum 
of the shock becomes softer as it expands and cools, its 
NIR light curve evolves on much shorter timescales than 
its redshifted bolometric luminosity. These events will be 
easily recognizable as transients and distinguished from 
protogalaxies. 

7. CONCLUSION 

Luminous Type Iln SNc will probe stellar populations 
at z = 10 - 20, redshifts that complement those at which 
normal-luminosity CC and PI SNc can be detected 
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Fig. 11— NIR light curves for all 5 SNe at z = 30. The dashed 
line at AB mag 32 again denotes the JWST NIRCam photometry 
limit. 

(z < 15 and z > 30, respectively). They will not 
trace the first generation of stars, those that form at 
z ^ 20 - 30, but they will be found during the rise 
of the Lyman- Werner background (10 < z < 20) and 
in the first galaxies at 2 ~ 10 - 15. The event rate 
of Type Iln SNe depends in part on what fraction 
of 20 - 40 Mq Pop III stars die as compact blue 
giants or shed a common envelop e in a binary, and 
may_ b e low (IWise fc Abe| ^ l2005t iWeinmann & : Lilryl 



2001 lO'Shea et al l 120051 iTornatore et al.l 12007 



Whalen et all I2008at iTrenti et all 120091; I Whalen et al 
20101: IGreif et al.ll20ld IMaio et al.ll2011t iHummel et al 



2012t Ijohnson et all 120131: IWise et all 1201 2D . For exam- 
ple. I Tanaka et alT i 20121) take this rate to be a few tenths 
of a percent of the total Type II SN rate. However, such 
estimates may be too conservative for higher redshifts 
because the Pop III IMF is known to be top heavy and 
because the evolution of 20 - 40 Mq primordial stars is 
still not fully understood. 

The other challenge to observing such explosions is 
that they are too dim to be detected beyond z > 7 in 
all-sky NIR surveys by WFIRST or WISH, whose wide 
fields of view would otherwise compensate for the low 
SN Iln event rate. This picture could change if the core 
collapse event is more energetic than the 2 .4 x 10 51 erg 
explos ion considered here. For example, iMoriva et al.l 
(|20H find that SN2006gy is best modeled by a 4 x 10 52 
erg hypcrnova explosion in a 15 Mq circumstellar shell. 
Such events, with or without shells, may be visible in 
all-sky NIR campaigns at much higher redshifts and will 
be the focus of a future paper. 

In our models we assume a uniform shell, but in reality 
it could be clumpy due to a variety of hydrodynamical 
instabilities. If so, the ejecta would light up the shell 
unevenly, with brighter emission emanating from denser 
clumps than from diffuse regions. It is not clear how the 
total luminosity of the collision would change in such cir- 
cumstances. However, this scenario would occur less of- 
ten in the primordial universe than today because there 
are no metals or dust to radiatively cool and fracture 
the shell. Clumping could still occur if the medium into 



which the shell is ejected is not uniform or if the outflow 
itself is collimatcd. The former is less likely at high red- 
shifts because the progenitor star usually drives all the 
gas from the halo in strong ionized flows, leaving behind 
diffuse uniform media in the vicinity of the star. 

Radio emission from these explosions may be visible 
at 21 cm by eVLA, eMerlin and the S quare Kilometer 
Array fffKAV iMeiksin fe Whaler] (|2012l ) found that syn- 
chrotron emission from CC SNe at z = 10 - 20 will be 
detected by SKA and that more energetic hypernovae at 
this epoch can be detected by existing facilities. Addi- 
tional calculations are necessary to determine if the col- 
lision of the ejecta with the shell enhances or quenches 
its radio emission. Type Iln SNe are not expected to im- 
print excess power on the CMB on small scales because 
unlike PI SNe, CC SNe are not sufficiently en ergetic to 
Comptonize large number s of CMB photons (|Oh et al.l 
l2^0ltlWhaTen et al.ll2008c^ . 

A small fraction of Pop III CC S Ne may proceed 
as g a mma ray bursts (GRBs; e.g., iBromm fc Loebl 
l2006t iWang et all 12012ft , driven eit her by the collapse 
of massive rapidly r otating stars (jSuwa fc Iokal l201lt 
INagakura et al.l I2012[) or binary mergers between 20 - 
50 Mq stars ( e.g.. IFrver fc Wooslevl 119981: IFrver et al.l 



[19991: IZhang fc Frverl 120011: IFrver et al.ll2007E This is 
reinforced by the fact that some Pop III stars have 
been found to form i n binaries in numerical simula- 
tions (ITurk et al.H2"00l . While x-rays from these events 
could be detected by Swift or next-generation mis- 
sions such as the Joint Astrophysics Nascent Universe 
Satellite ( JANUSM&^&ios fc Reesl l20Ul IRomindl200l 
iBurrows et al.l l2010l). it is m ore likely that their af- 
terglows (|Whalen et al.l l2008bf ) will be detected in all- 
sky radio surveys by the eVLA, eMERLIN and SKA 
(|de Souza et al.l [20 111) due to their low event rates. If 
these events occur in dense circumstellar shells like those 
in our models, the s hells may imprint d istinct features on 
the afterglows (e.g., iMesler et al.ll2012l) . We arc now de- 
termining the afterglow signatures of Pop III GRBs in a 
varie ty of circumburst environments (Mesler et al.ll2012l 
|20T1 . 

Pop III Type II SNe will completely outshine the 
primeval galaxies in which they occur because they have 
comparatively few stars. These events may reveal the 
existence of such galaxies when they might not other- 
wise have been detected by JWST or future 30 m class 
telescopes such as the Giant Magellan Telescope or the 
Thirty-Meter Telescope. Together with CC and PI SNe, 
Pop III shell-collision SNe will trace star formation rates 
and chemical enrichment in nascent galaxies. These an- 
cient explosions will soon open a new window on the 
high-rcdshift universe. 
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